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SMART ANTENNA SYSTEM WITH IMPROVED LOCALIZATION OF 

POLARIZED SOURCES 

5 

FIELD AND BACKGROUND OF THE INVENTION 
The present invention relates to a smart antenna method and apparatus and, 
more particularly, but not exclusively to the use of polarization as an additional 
parameter for determining the angle of arrival of impinging electromagnetic radiation. 

10 Vector sensors enable estimation of the angle of arrival and polarization of 

impinging electromagnetic waves with arbitrary polarization. In the last decade, many 
array processing techniques for source localization and polarization estimation using 
vector sensors have been developed. A. Nehorai and E. Paldi, "Vector-sensor array 
processing for electromagnetic source localization," IEEE Trans, on Signal 

15 Processing, vol. 42, pp. 376-398, Feb. 1994, the contents of which are hereby 
incorporated by reference, developed the Cramer-Rao bound (CRB) and the vector 
cross-product direction-of-arrival (DOA) estimator. Polarimetric modeling using 
vector sensors is performed in B. Hochwald and A Nehorai, "Polarimetric modeling 
and parameter estimation with application to remote sensing" IEEE Trans. On Signal 

20 Processing, Vol. 43 pp. 1923-1935, Aug 1995 the contents of which are hereby 
incorporated by reference. 

Identifiability and uniqueness issues associated with vector sensors are 
analyzed in G. F. Hatke, 27 th Asilomar Conf. 1993, pp 1365-1369, K. C. Ho, K.C. 
Tan, W. Ser, Signal Processing Vol. 47 pp 41-54 Nov. 1995 , B. Hochwald and A 

25 Nehorai, IEEE Trans on Signal Processing, Vol. 44 No. 1 pp 83-95 Jan 1996, K. C. 
Tan, K.C. Ho, A. Nehorai IEEE Trans. Signal Processing, Vol. 44 pp. 3099-3107, 
Dec. 1996. The contents of each of the above documents are hereby incorporated by 
reference. 

Eigenstructure-based techniques, such as ESPRIT and multiple signal 
30 classification (MUSIC) for source localization using vector sensors have been 
extensively investigated. J. Li, IEEE Trans. Antenna Propagation. Vol. 41. pp 379- 
387 March 1993, the contents of which are hereby incorporated by reference, applied 
the ESPRIT algorithm to a vector sensor array. ESPRIT-based direction finding 
algorithms using vector sensors have been further investigated in several papers 
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including K. T. Wong and M. Zoltowski, IEEE Trans. Antenna Propagation. Vol. 48 
pp 671-681, May 2000, the contents of which are hereby incorporated by reference. 
MUSIC-based algorithms for the same problem have been applied in K. T. Wong and 
M. Zoltowski, IEEE Trans. Antenna Propagation. Vol. 48 pp 1235-1245 and 2205- 

5 2210, August 2000. These techniques yield high-resolution and asymptotically 
efficient estimates in case of uncorrelated or partially correlated signals. However, 
since these techniques assume a non-singular signal correlation matrix, they encounter 
difficulties in cases of fully correlated signals such as signals from a single source in a 
multipath scenario. The latter is particularly common in cellular telephony in the 

10 urban environment. 

In order to decorrelate signals in a data covariarice matrix, Evans et al. in Proc. 
1 st ASSP Workshop spectral Estimation, Hamilton Ontario Canada, 1981 ppl34-139, 
the contents of which are hereby incorporated by reference, proposed a preprocessing 
technique referred to as spatial smoothing. Several later authors investigated the 

15 spatial smoothing method, in combination with the method of forward-backward 
averaging. The drawback of the above combination is the reduction of the effective 
array aperture length, resulting in lower resolution and accuracy. An alternative 
spatial averaging method is redundancy averaging. It has been shown that a 
preprocessing method based on redundancy averaging induces bias in the DO A 

20 estimates. 

A Maximum Likelihood (ML) approach for diversely polarized source 
localization was proposed in Wax and Ziskind, IEEE Transactions Antenna 
Propagation vol 38, pp. 111-1114, July 90. The maximum likelihood method uses the 
simulated annealing algorithm to efficiently search over the signal DOA's and 
25 polarization vectors. 

The above approach however still does not overcome the difficulties 
introduced by having correlated signals, as would be expected when seeing multi- path 
versions of the same signal. 

There is thus a widely recognized need for, and it would be highly 
30 advantageous to have, a decorrelation method usable with a smart antenna, which is 
devoid of the above limitations and is therefore useful for cellular telephony in urban 
areas and for other cases where multi-path is a significant problem. 
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SUMMARY OF THE INVENTION 

According to one aspect of the present invention there is provided a quadrature 
polarized antenna element comprising a plurality of electric dipoles arranged at a 
predetermined angle with respect to one another and a plurality of magnetic dipoles 
arranged at said predetermined angle with respect to one another, and wherein each 
magnetic dipole substantially shares a common location with a respective one of said 
electric dipoles. Preferably, said predetermined angle is substantially a right angle. 
Furthermore, said plurality of dipoles are arranged to obtain electro-magnetic source 
information from throughout the azimuth plane. 

Preferably, said electromagnetic source information is such as to allow 
extraction of at least one of a group comprising signal localization, signal polarization, 
and data content. 

The element is preferably arranged to detect signals from substantially 
throughout the azimuth plane, said electric dipoles and said magnetic dipoles being 
arranged to complement each other in respect of polarization information in said 
signals, thereby to enable said element to obtain substantially all polarization 
information within said azimuth plane. 

Preferably, said element comprises four dipoles, of which two of said dipoles 
are electric dipoles, and two of said 'dipoles are magnetic dipoles. In a preferred 
embodiment, the two electric dipoles are mutually orthogonal and said magnetic 
dipoles are mutually orthogonal, and each of said magnetic dipoles is co-directed with 
one of said electric dipoles. 

According to a second aspect of the present invention, there is provided a 
quadrature polarized antenna array, comprising a plurality of antenna elements, each 
element comprising a plurality of electric dipoles arranged at a predetermined angle 
with respect to one another and a plurality of magnetic dipoles arranged at said 
predetermined angle with respect to one another, and wherein each magnetic dipole 
substantially shares a common location with a respective one of said electric dipoles. 
The properties of the individual elements are as described above, and the arrangement 
of the elements into an array provides additional properties to the array such as the 
ability to resolve originating sources from multipath variations of the signal. For an 
array having n elements, the array is able to resolve n-1 signals. 



In one embodiment, the array comprises an electric switch or switching 
element for switching between dipoles or a switch for switching between elements to 
gather data over said array. 

In one embodiment, the array has connected thereto a signal preprocessor for 
preprocessing signals from said antenna for obtaining spatial spectrum information for 
signal source location, the preprocessor comprising: 

a sensor autocorrelator configured for forming signal autocorrelation matrices 
for each sensor type, and 

a smoother configured for smoothing said autocorrelation matrices, thereby to 
form at least one covariance matrix comprising spatial spectrum information. 

Preferably, the antenna array has connected thereto a source locator configured 
for using said sample covariance matrix in an eigenstructure-based signal source 
localization technique. 

Preferably, said smoother is further configured to apply forward backward 
smoothing to said sample covariance matrix, thereby to increase a maximum number 
of signal sources that can be localized. 

The array preferably has a steering vector of antenna specific parameters for 
use together with said sample covariance matrix in said source locator. 

In one embodiment, the array has connected thereto a signal preprocessor for 
preprocessing signals from said antenna for obtaining spatial spectrum information for 
signal source location, the preprocessor comprising: 

an autocorrelator, connected after said input, for forming signal autocorrelation 
matrices for each element, and a smoother, connected after said autocorrelator, for 
smoothing said autocorrelation matrices, thereby to form a sample covariance matrix 
suitable for use in a Eigenstructure-based estimator for estimating source localization. 

According to a third preferred embodiment of the present invention there is 
provided a method for preprocessing incoming signals obtained using a plurality of 
different sensor types, the signals including coherent signals, the preprocessing being 
for source localization, the method comprising obtaining angle of arrival and 
polarization information of incoming signals from each of said different sensor types, 
forming signal autocorrelation matrices for each sensor type, and smoothing said 
autocorrelation matrices, to form therefrom at least one covariance matrix suitable for 
use in eigenstructure-based signal source localization techniques. 



5 

The method may further comprise applying forward backward averaging to 
said sample covariance matrix, thereby to increase a maximum number of signal 
sources that can be localized. 

The method may further comprise obtaining, at least once for each antenna, a 
steering vector for use together with said sample covariance matrix in said 
eigenstructure-based signal localization techniques. 

Sensing or obtaining in the method is from four sensor types, Preferably, the 
sensor types are two respectively orthogonal electrical dipoles and two respectively 
orthogonal magnetic dipoles. 

Preferably, the four sensor types are all arranged for sensing in a single plane. 
That is to say they are arranged to give full coverage of that plane, but of course they 
are still able to sense signals coming from above or below the plane. 

The method may comprise using source localization information obtained from 
the data of said covariance matrix as an input to a beam director to provide a directed 
beam to a respective source. 

Preferably, said incoming signal is a noise signal, and the method may 
comprise using source localization information obtained from said covariance matrix 
as an input to a beam director to provide a null of a directed beam to a respective 
source of said noise interference signal. 

According to a fourth aspect of the present invention there is provided a 
method for processing incoming signals obtained using a plurality of different sensor 
types, the signals including coherent signals, the preprocessing being for source 
localization, the method comprising obtaining angle of arrival and polarization 
information of incoming signals, forming signal autocorrelation matrices over an array 
of said sensors, and forming a sample covariance matrix from said signal 
autocorrelation matrices, said sample covariance matrix being suitable for use in a 
maximum likelihood estimator for estimating source localization. s 

Preferably, said maximum likelihood estimator is 



(^) = argmaxA taax 

8,$ 



wherein: 

6 - vector of the elevation angles of the sources, 



6 

<j> - vector of the azimuth angles of the sources, 

/L inax - maximum generalized eigenvalue of the matrix pair (*F 2 (#,^),*F,(#,^)) 
F T (&,0) - matrix whose columns denote spatial transfer functions for both polarization 
components of the sources, 

H - matrix hermitian operation (complex conjugate and transpose) 
R„ - noise and interference co variance matrix 

R^- sample covariance matrix 

The method may use source localization information obtained from said 
covariance matrix as an input to a beam director to provide a directed beam to a 
respective source. 

Typically, the incoming signals includes noise, for which the method further 
comprises using source localization information obtained from said covariance matrix 
as an input to a beam director to provide a null of a directed beam to a respective 
source of detected noise. 

According to a fifth aspect of the present invention there is provided apparatus 
for preprocessing incoming signals obtained using a plurality of different sensor types, 
the signals including coherent signals, the preprocessing being for source localization, 
the apparatus comprising an input for obtaining angle of arrival and polarization 
information of incoming signals from each of said different sensor types, a sensor 
autocorrelator configured for forming signal autocorrelation matrices for each sensor 
type, and a smoother, configured for smoothing said autocorrelation matrices , thereby 
to form at least one covariance matrix suitable for use in eigenstructure-based signal 
source localization techniques. 

Preferably, said smoother is further configured to apply forward backward 
averaging to said sample covariance matrix, thereby to increase a maximum number of 
signal sources that can be localized. 

The apparatus may make use of a steering vector for use together with said 
sample covariance matrix in said eigenstructure-based signal localization techniques. 

The apparatus preferably makes use of the four sensor types discussed above. 

The apparatus may comprise a covariance unit, also connected after said 
autocorrelator, for forming a covariance matrix from which a sample covariance 
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matrix \can be extracted, said sample covariance matrix being suitable for use in a 
maximum likelihood estimator for estimating source localization. 

The apparatus may comprise a switching unit for switching between said 
smoother and said cross-correlator in accordance with a total number of apparent 
signal sources giving rise to said incoming signals. 

According to a sixth aspect of the present invention, there is provided 
apparatus for processing incoming signals obtained using a plurality of different sensor 
types, the signals including coherent signals, the processing being for source 
localization, the apparatus comprising 

an input for obtaining angle of arrival and polarization information of incoming 
signals, an autocorrelator, connected after said input, for forming signal 
autocorrelation matrices for each element, and a covariance unit connected after said 
autocorrelator for forming a covariance matrix, from which a sample covariance 
matrix is extractable from said autocorrelation matrices, said sample covariance matrix 
being suitable for use in a maximum likelihood estimator for estimating source 
localization. 

Preferably the maximum likelihood estimator is defined as 



wherein: 

6 - vector of the elevation angles of the sources, 
(f> - vector of the azimuth angles of the sources, 

/l max - maximum generalized eigenvalue of the matrix pair (?¥ 2 (0,0\*¥ l (0 9 (fi)) 

¥ T {0,<f) - matrix whose columns denote spatial transfer functions for both polarization 

components of the sources, (defined in (20)) 

H - matrix hermitian operation (complex conjugate and transpose) 

R„ - noise and interference covariance matrix 

R - sample covariance matrix 



Unless otherwise defined, all technical and scientific terms used herein have 
the same meaning as commonly understood by one of ordinary skill in the art to which 
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this invention belongs. The materials, methods, and examples provided herein are 
illustrative only and not intended to be limiting. 

Implementation of the preprocessor and preprocessing method and system of 
the present invention involves performing or completing selected tasks or steps 
5 automatically using software or hardware or a combination thereof. Moreover, 
according to actual instrumentation and equipment of preferred embodiments of the 
method and system of the present invention, selected steps may be implemented in 
hardware or software on any operating system or any firmware or a combination 
thereof. For example, as hardware, selected steps of the invention could be 
10 implemented as a chip or a circuit. As software, selected steps of the invention could 
be implemented as a plurality of software instructions being executed by a computer 
using any suitable operating system. In any case, selected steps of the method and 
system of the invention could be described as being performed by a data processor, 
such as a computing platform for executing a plurality of instructions. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to 
the accompanying drawings. With specific reference now to the drawings in detail, it 
is stressed that the particulars shown are by way of example and for purposes of 

20 illustrative discussion of the preferred embodiments of the present invention only, and 
are presented in the cause of providing what is believed to be the most useful and 
readily understood description of the principles and conceptual aspects of the 
invention. In this regard, no attempt is made to show structural details of the invention 
in more detail than is necessary for a fundamental understanding of the invention, the 

25 description taken with the drawings making apparent to those skilled in the art how the 
several forms of the invention may be embodied in practice. 
In the drawings: 

FIG. 1 is a schematic diagram illustrating two cellular base stations and 
showing how source location according to preferred embodiments of the present 
30 invention lead to more efficient use of the spectrum in the case of a single base station 
tracking two moving signal source/targets; 

FIGS. 2a - 2d are simplified diagrams illustrating coverage patterns for mobile 
base stations in ideal and real environments; 
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FIG. 3 is a simplified schematic illustration of an electrical dipole antenna; 
FIG. 4 is a simplified schematic illustration of two electrical dipoles arranged 
orthogonally; 

FIG. 5 is a simplified schematic illustration of a magnetic dipole; 
5 FIG. 6 is a simplified schematic illustration of two magnetic dipoles arranged 

orthogonally; 

FIG. 7 is a simplified schematic illustration of an antenna unit comprising two 
orthogonally arranged electrical dipoles and two orthogonally arranged magnetic 
dipoles placed together according to a first preferred embodiment of the present 
10 invention; 

FIG. 8 is a simplified diagram showing a first preferred embodiment operative 
in accordance with the present invention, of a method for obtaining spatial information 
from an array of elements such as the element of Fig. 7; 

FIG. 9 is a simplified flow chart illustrating an alternative preferred 
15 embodiment, operative in accordance with the present invention, of preprocessing to 
obtain spatial information from an array of elements; 

FIG. 10 is a simplified schematic diagram illustrating the advantages of 
accurate source location according to embodiments of the present invention, in a two 
source two base-station system with multiple path errors; 
20 FIG. 1 1 is a simplified diagram illustrating the antenna array unit of Fig. 7 laid 

out on a three dimensional grid; 

FIG. 12 is a simplified diagram illustrating switching between the antenna 
elements; 

FIG. 13 is a graph that compares results based on a batch of 100 data samples 
25 (snapshots), for different decorrelation methods including those in accordance with the 
present invention; 

Fig. 14, is a further graph showing simulation results, and shows RSME versus 
angle rj 2 at SNR = 20dB ; 

Fig. 15 is a further graph showing simulation results, in particular the case of 4 
30 coherent sources arriving from directions -30°,-15°,0 0 and 20° having different 
polarizations; 
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Fig. 16, is a further graph showing simulation results, in this case showing 
results using four additional coherent sources from directions -60°, -50° ,30° and 

45° and giving a total of eight coherent sources which are successfully identified as 
being coherent by embodiments of the present invention but not by the prior art. 

5 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present embodiments describe a smart antenna and associated methods for 
source localization using array vector sensors. In one embodiment a maximum- 
likelihood estimator provides source location for fully correlated sources. In a second 

10 embodiment, a preprocessing method overcomes the problem represented by 
singularities in the signal correlation matrix. The second embodiment is based on a 
procedure referred to herein as vector sensor smoothing (VSS), which provides spatial 
spectral information for the use of spectral based algorithms such as the eigenstructure 
based MUSIC, MVDR and ESPRIT methods, for DOA estimation in scenarios 

15 involving fully correlated signals. Once the sources have been located, it is possible 
to communicate using directed beams and thereby to make more efficient use of the 
spectrum. The present methods allow dynamic real time source location so that 
directed beam technology can be used for moving sources. 

Also disclosed is a quadrature polarized antenna array optimized for obtaining 

20 polarization information in the azimuth plane. Use of any of the above methods 
together with the quadrature polarized array provides improved location performance 
compared to use of the respective method alone. 

The principles and operation of a smart antenna according to the present 
invention may be better understood with reference to the drawings and accompanying 

25 descriptions. 

Before explaining embodiments of the invention in detail, it is to be understood 
that the invention is not limited in its application to the details of construction and the 
arrangement of the components set forth in the following description or illustrated in 
the drawings. The invention is capable of other embodiments or of being practiced or 
30 carried out in various ways. Also, it is to be understood that the phraseology and 
terminology employed herein is for the purpose of description and should not be 
regarded as limiting. 
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Referring now to the drawings, Fig. 1 is a schematic diagram illustrating two 
cellular base stations and illustrating an advantage of the present embodiments. A first 
mobile base station 10 is in communication with a range of cellular telephones 12, 14 
and 16. The base station 10 is unable to determine with great accuracy the direction 
5 from which the signals come due for example to the multipath effects of signals 
bouncing off buildings. Thus, the signals are not strongly directed. Reference 
numeral 18 indicates the space / power distribution of the various signals. As is 
apparent from the illustration, the direction of the secondary node of the signal to 
source 12 coincides with the direction of the main node to source 14. Such prevents 

10 sources 12 and 14 from using the same or otherwise interfering frequencies. It is not 
possible to use a narrower beam for either of the sources because the directions of the 
sources are not known accurately. 

Base station 20 is able, in accordance with a preferred embodiment of the 
present invention, to determine with accuracy the direction of sources 22 and 24, as 

15 will be explained in detail below. Not only is it able to locate static sources but, as the 
location can be achieved in real time, even dynamic sources 22 and 24 can be tracked. 
Source 24 is a vehicle, which it is possible to track with a narrow beam, according to a 
preferred embodiment of the present invention. 

Reference is now made to Fig. 2, which is a simplified diagram illustrating 

20 specific problems of cellular telephone base stations in the urban environment. Fig. 2a 
illustrates the coverage area of a base station over smooth uninterrupted terrain. Fig. 
2b shows the coverage area of a directional antenna over smooth interrupted terrain. 
Fig. 2c shows the coverage area of a base station over urban terrain in which buildings 
provide obstacles to the free passage of signals. Fig. 2d shows the coverage area of a 

25 directional antenna over suburban terrain, that is to say with a lower density of 
building than that of Fig. 2c. Buildings in the urban and suburban environment cause 
signals to be blocked or reflected, and thus make signal source location difficult. The 
cellular base station may see the same signal approaching from several different 
directions, a phenomenon known as multipath error, or may not see the signal arrive at 

30 all. Current algorithms for source location tend to fail beyond two multipath errors. 
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The Antenna 

Figs. 3-7 show sensing antenna dipoles, and how they may be combined into a 
quadrature polarized antenna according to a first embodiment of the present invention. 
Reference is now made to Fig. 3, which illustrates a typical electrical dipole. 
5 The dipole 30 comprises two laterally extending poles 32 and 34. In Fig. 4, two such 
dipoles 30 are placed together orthogonally. 

Reference is now made to Fig. 5, which shows a typical magnetic dipole. The 
dipole 40 is formed by a single wire describing a rectangle. 

In Fig. 6, two magnetic dipoles are shown placed together orthogonally. 
10 Reference is now made to Fig. 7, which is a simplified diagram illustrating an 

antenna element according to a preferred embodiment of the present invention. Two 
electric dipoles 30 are placed together orthogonally as illustrated in Fig. 4. In addition, 
two magnetic dipoles 40 are also placed together orthogonally as shown in Fig. 6, 
together with the electrical dipoles to form an antenna element having four mutually 
15 orthogonal sensing dipoles. All four dipoles sense in the same plane. If the antenna is 
oriented as shown in the illustration then sensing is maximized in the azimuth or 
horizontal plane, and the sensors allow polarization information within the azimuth 
plane to be obtained. As may be appreciated, when considering cellular 
communication, most of such communication is limited to the azimuth plane. Cellular 
20 signals are polarized and thus polarization may be used as an additional parameter, as 
will be explained below, to eliminate multipath errors. Deviations from the azimuth 
plane can also be sensed although substantially vertical signals may be lost. 

Vector Sensor Smoothing (VSS) 

25 Reference is now made to Fig. 8, which is a simplified diagram 

showing a first preferred embodiment of a method for obtaining spatial information 
from an array of elements such as the element of Fig. 7. In a stage SI, information is 
obtained from the various sensors and elements of the array. The data in each sensor 
is not exactly the same. That is to say for multipath versions from the same signal 

30 source the data carried by the signal remains the same but the specific parameter 
picked up by each individual dipole is different. If the dipoles are arranged to be 
orthogonal as in Fig. 7, and if the signal is polarized, then the amount of information is 
maximal. Then, in a stage S2, autocorrelation matrices are formulated for each of the 
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sensor types, that is dipole configurations, in the array. That is to say, one 
autocorrelation matrix is formulated for electrical dipoles in a first sense, a second 
autocorrelation matrix is formulated for electrical dipoles in the orthogonal sense. A 
third autocorrelation matrix is formulated for magnetic dipoles in a first sense and a 
5 fourth autocorrelation matrix is formulated for magnetic dipoles in the orthogonal 
sense. 

In a stage S3, an average is taken of the autocorrelation matrices, which is to 
say that the autocorrelation matrices are summed and then divided by the total number 
of matrices. The result is an extended covariance matrix and is referred to 
10 hereinbelow as R VSS5 where vss refers to vector sensor smoothing. 

An optional stage S4 follows of applying forward backward averaging to the 
matrix R vss to form R VS s-fb- Forward backward averaging is a method, known to the 
skilled person, for decorrelating signals, and reference was made thereto hereinabove. 
More particularly, forward backward: averaging is a process applied to the 
15 autocorrelation matrix using the complex conjugate of the autocorrelation matrix 
multiplied by anti-diagonal matrixes of both sides thereof. Mathematically forward 
backward averaging is carried out as follows: 

20 where: 



0—0-1 
! 1 0 

1 0—0 

The use of forward backward averaging is advantageous in that the maximum 
25 number of signal sources that can be localized is doubled. However, the forward 
backward averaging method assumes a symmetric array, far field approximation and 
unequal signal phases at the center of the array. Thus it is not always applicable. In 
general, for mobile telephone sources these assumptions, in particular that of the 
symmetric array, are true. If the assumptions are not true then it is advisable to use 
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vector sensor smoothing without forward backward averaging. Other methods of 
decorrelating signals may also be considered. 

In stage S5, a steering vector q, as defined in (14) below, is used. Although 
shown sequentially, stage S5 is independent of the immediately preceding stages and 
5 thus may be carried out in parallel to save processing time. 

In stage S6, R vss -fb if available from stage S4, or otherwise R VSS5 is passed on for 
the use of spectral based algorithms such as MUSIC and MVDR. The steering vector 
q is also passed on. Theoretically a steering vector could be generated from any of the 
autocorrelation matrices, that is to say for each of the sensor elements. Preferably, 
10 however, just a single steering vector q is passed on. Since all of the sensors obtain 
the same signal there is no significant difference between the steering vectors obtained 
from the different sensors. 

The procedure described in stages 1-5 serve as a preprocessing stage for the 
spectral based algorithms, which then complete the task of characterizing multipath 
15 versions of the same signal. 

Maximum Likelihood Estimation 
Reference is now made to Fig. 9, which is a simplified flow chart illustrating 
an alternative method of obtaining spatial information from an array of elements of 

20 the kind described above. In Fig. 9, a first stage S10 obtains information from the 

dipoles or sensors of an antenna array element of the kind shown in Fig. 7. In a stage 
SI 1, the information is autocorrelated over the individual elements, just as was done 
in the method of Fig. 8, to form autocorrelation matrices R n j. In a stage SI 2, a 
sample co variant matrix A Ry is formulated. In a stage S 13 the sample co variant 

25 matrix is supplied as an estimate to a maximum likelihood estimator. In a stage S14 
the maximum likelihood estimator is used to obtain the direction of the source within 
the azimuth plane. 

Comparison Between Methods 
30 Maximum likelihood (ML) estimation achieves a greater accuracy than the 

VSS method of Fig. 8 and the Eigenstructure-based techniques, and approaches the 
Cramer Rao bound for maximum theoretical resolution. However ML is 
computationally expensive, especially for large numbers of incoming signals. For 
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small numbers of incoming signals, including one, however, the computational cost is 
not that great and can be justified by the increased accuracy. Thus a preferred 
embodiment of a direction finder according to the present invention incorporates both 
a VSS estimator according to Fig. 8 and a maximum likelihood estimator according to 
5 Fig. 9 and chooses a suitable threshold of signal paths at which to swap between the 
two methods. 

In any of the above embodiments it is possible to provide dedicated hardware 
for gathering data from individual sensors and individual elements in the array, and 
building the autocorrelation matrices therefrom. In an alternative embodiment it is 
10 possible to take) data from the different sensors in turn or from different elements in 
turn using high-speed switching techniques. In a particularly preferred embodiment a 
single data gathering unit is used for each of the four sensors in each element and 
switching is provided between the sensors. Thus the antenna saves on hardware. 

15 Cellular Telephony Using the Present Embodiments 

Reference is now made to Fig. 10, which is a schematic diagram illustrating the use of 
source location in the cellular telephony environment. As the base station is able to 
localize the source, the same channel can be reused in different directions, even where 
the directions are fairly close, and can certainly be reused by neighboring base 

20 stations. In Fig. 10, source 100 transmits a signal to base station 102. The signal 
reaches the base station via five different paths. Using the methods described above, 
the base station 102 is able to identify that a single source is responsible, locate the 
source and communicate with it using a relatively accurately directed beam. Source 
104 is then able to use exactly the same channel to communicate with neighboring 

25 base station 106. The signal also manages to reach base station 102. Base station 102 
does not mistakenly identify the signal with source 100 since the polarization is 
different. Thus no channel interference is created. 

Such an antenna is also able to cancel a noise source by directing a null of a 
given channel towards an identified source of noise. 

30 In the following, a mathematical description of the embodiments is provided. 

A measurement model using a vector sensor array is discussed. Then an ML estimator 
for DOA and polarization vectors is derived. The derivation is followed by the vector 
sensor smoothing (VSS) method as a preprocessing stage for eigen-structure based 
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source localization. The performance of the proposed algorithms is evaluated via 
computer simulation and described. 

Signal Characterization 
We begin by considering the polarization of an incoming signal. Assume the 
5 incoming signal to be a planar wave incident from the direction {9*<f>\ where 6 is the 
elevation angle and <j> is the azimuth angle. The complex envelope of the transverse 
incident electric field can be expressed in spherical coordinate system by 

(1) E = Ej) + E e Q 

10 

where and E e are its horizontal and vertical components. For a given 

polarization, the electric field components at the sensor can be described by 

15 (2) E^(t) = Acosys(t) 

and 



(3) E e (t) = Asiny e jn s(t) 

20 

where A is the signal complex amplitude. The parameters y and rj determine the 
polarization and are related to the amplitude and phase of the horizontal and vertical 
components of the received signal. s(t) is a base-band equivalent of the transmitted 
narrow band signal. 

25 

The electric field in (1) can be expressed in x,y 9 z Cartesian coordinate system by 

E = ^E e cos0 cos^ -i^sin^ )x 

(4) . +(i?0cos# sin^ +E <p cos<f> )jp 

-(i^sin^? )f 
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By substitution of (2) and (3) into (4) one obtains 



E(t) = A^sin ye jn cos 9 cos <f> — cos y sin <f> )x 
+ {sin ye Jn cos 6 sin <f> + cos y cos (p )y 
-{siny e Jt1 sin0)z)p(t) 



(5) 



5 Any plane wave can be characterized by the angles of arrival 0,<f>, the complex 
amplitude A , and the polarization parameters y and 77 . 

The polarization parameters can be expressed by vector notation as 



10 



pi 




sinye Jn 


(6) 


p*. 




cosy 
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P = 

and E(t) can be rewritten in the form: 

E(t) = As(t){p d [xcos0cos<f> + ycos6sin(/>-zsin0 ]+pj[-xsin0 + y cos (7) 
Similarly, the incident magnetic field is given by 



£(0 = 



in which Z 0 is the characteristic impedance of the medium. 



(8) 
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Array Vector Sensor Spatial Response to a Polarized Signal 

1 . Sensor Response 
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We now consider a general type of monopole or dipole sensor, which may be either 
electric or magnetic and linearly polarized. Let V x9 V y9 V z denote the sensor response 

to x, y, z, components of an arbitrary incident electric or magnetic field. Accordingly, 
the total signal at the sensors output terminals could be described by 

gAO^v) = a E (0^)p = [a e {0 f if) a,{0,<f>)] 

and 

gn (#> 0, P) = "h ( 0 > <f>)P = k *) ~ a o & *)] 
in which the subscripts E and //stand for the electric and magnetic fields and 

a e {0,</>) = V x cosO cos<p+V y cosO sin<f>-V z sinO 

(10) . 

a <t> (&> <f>) = ~ V x Sin <f>+ V y COS <f> 



Pe 



(9) 




2. Vector Sensor Response 
Reference is now made to Fig. 11, which shows a multiple orthogonal 
component sensor for probing both magnetic and electric elements in the three 
dimensional plane. 

Consider a vector sensor containing 3 electric and 3 magnetic orthogonal 
sensors, located at (x 09 y 0 ,z Q ) and azimuthally rotated by an angle 8 , as depicted in 

Fig. 1 1 . The coefficients V x , V and V z for these sensors are 



electric sensors: 



magnetic sensors: 



V xl =sm(5) ^=-cos(J) V 2l =0 
V x2 = cos (S) V y2 = sm(S) V z2 = 0 



r y3 = Q 



^3=1 



V xA = wcos{S) V y4 = wsin(S) V zA = 0 
V x5 = wsin(S) V y5 = -wcos(5) V l5 = 0 



(ii) 
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where w denotes the ratio between the induced voltage in an electrical sensor to the 
corresponding induced voltage in a magnetic sensor. 

The spatial response in matrix notation of the vector sensor, shown in Fig. 1 1 , can be 
expressed by 



g 0 (^^P) = 



Pe 



(12) 



in which a^ k (0 s ^) 9 and a 0k (0,<p) are defined in Eq. (10) for the £th sensor with the 
corresponding V xk 9 V yk and V zk . The polarization p is determined by two real 
parameters, y and r\: 



f • in 
cos/ 



as in equation (6) above 



In general, a vector sensor may contain part of the six elements described hereinabove 
and therefore the corresponding spatial response vector size is given by 1<L<6. 

3. Array Vector Sensor Response 
For the general case of a 3D array with N vector sensors, the spatial response in 
matrix notation of the array vector sensors is expressed by 

gfo ^ P) = q(0, ft ® go (ft & P) - [q(g» 4>) g A(0, gj p (13) 
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where ® denotes the Kronecker product The size of the vector q(#,^) is jVxland 
its elements represent the phase delay associated with each vector sensor in the array 
due to its relative location for an incident plane wave from the direction (0,<f>): 

5 q(M = 

and k Q is the wavenumber in the medium. 

The Measurement Model 

10 Consider the scenario of M signals, s m (t), impinging on the array from directions 
(0 m ,<f> m ) and polarization vector p m , where m =T, . . . , M . The spatial response of the 
array to the mth signal is denoted by g{0 m ,& m ,p m ) . Assuming K independent samples 
of the measurement vector, y(f ) ,the data model is given by 

15 y(k) = f i F(0 m ,0 m ,p m )s mk +n k (15) 

where K is the number of independent samples collected by the array and n k 
represents the kth sample of the additive noise and interference vector. 
The measurement and noise vectors, y*and are each of size LN , the matrix 

20 F(0 W ,^ OT ) is of size LNx2 whose columns denote the spatial transfer functions for 
both polarization components of the mth signal, and p m is a complex vector of size 2 
describing the corresponding signal polarization state. 

We assume that the noise vector, {n^}f =l , is an Li.d. sequence with zero-mean, 
complex Gaussian distribution, n k - CA^(0,R n ), and independent of the signals. The 
25 signals snapshots, s mk , are assumed to be unknown deterministic. In case of coherent 
or correlated signals, they can be decomposed as 



g J k o 1*1 sin 0 cos # + yi sin e sin <f>+Z\COS0 \ 
Jk 0 [x 2 sin 0 cos <p+y 2 sin 0 sin <f>+z 2 cos0\ 



(14) 



Jk 0 [xtf sin 0cos$+y N sin 0 sin#'+z N cos0 ] 
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s mk = ji m s k9 w = 2,...,M, fc = l,...,AT, where // m denotes the relative amplitude 
and phase of the mth signal. Thus, Eq. (15) can be rewritten in the form: 

M 

y k =Tv{O m ,<f> m h m » m s k +" k , k = \,...,K. (18) 

m=\ V v ' 

5 

By denoting g m =p m // m , the unknown parameters space is reduced to the source 

directions 0 - (0 } , . . . , 6 M ) r , ^ = {(f> { , . . . , ^ ) r , the modified signal polarization vector, 

g = {g? ,.~ 9 Glt) T , and the signal s =(si..., s k ) T . Eq. (18) can be rewritten in short form 
notation as 

10 

y* = F t(&i</>)$ s k + n k> k = l,...,K 

where 

15 

(19) g = [g[ gl ... g T M f 

and 

20 (20)F r (ft^)=[F(^) F(0 2 ,<* 2 ) ... F(0 M J M )] 

The problem of localizing the signal source is to estimate the directions of arrival, 
(0 m9 $ m ), while the signal vector s = (s l9 ...,s K ) T and the vector of the modified 
25 signal polarizations of the M arrivals q are unknown complex vector parameters. 
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Algorithm 



The Maximum Likelihood Estimator 

In this section, the ML estimator is derived for the problem stated in the previous 
section. Under mild regularity conditions, the ML estimator asymptotically achieves 
5 the Cramer-Rao Lower Bound (CRLB), as the number of snapshots goes to infinity. 



The ML estimator of the source location is given by 

0,^ = argmaxk (21) 

where L y ^) is the localization function, defined as 

L y (0j)=maxL y (0 9 f 9 g\s) (22) 



10 



25 



15 andZ^(6^,^,s) = log/(y 15 ...y^/#,0,£\s) is the conditional probability density 
function (pdf) of the measurements y^^y^ given the unknown parameters 0,<f>,g 
and s . 

Considering the assumptions stated in the previous section, y l9 ...,y K is an i.i.d. 

f \ 



20 sequence with y^ ~ CN 



Accordingly, 

The ML estimator of the signal vector, s, can be obtained by equating the 
corresponding derivative of the log-likelihood function to zero, which results in 

s k =k"K:%y$"R- n l y k (24) 
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Substitution of (24) into (23) yields 



l>, — * io g (^|R n |) - X y"Ky„ + * q ^ p -./ • < 26) 



5 By denoting: 

and using the definition ^ = F T (0,<fi)g , we obtain a simplified form of Eq. (26): 

10 

Ly(y „..., y ^,^)-C + K^^ (27) 

Maximization of the log-likelihood function in Eq. (27) with respect to the modified 
polarization vector, <; , is equivalent to finding the following maximum generalized 
15 eigenvector and corresponding eigenvalue: 

{y 2 (0,<j>)u = hv x {e^)u} (28) 

Thus 

^(^^^s) = C + iTA_{ M / 2 (^^), M / l (^^)} (29) 

20 

(29.1) $ = u^(v 2 (0j),v l 0J)) 

where ^ max {v 25 Vi} anc * ^maxi^'Vi} stand for the maximum generalized eigenvalue 
and the corresponding eigenvector of the matrix pencil (\|/ 2 (^,^),\j/ 1 (^ ? ^)) 3 
respectively. 

25 

Finally, the ML estimator of the signal DOA's is given by 
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0, j>) = arg max A max (\j/ 2 (<9, <*>), (0, #)) 



(30) 



and the estimate of the polarization vectors is given by the corresponding generalized 
eigenvector, 



The ML estimator is known to be asymptotically optimal, however it involves a 
2M dimensional search procedure for estimating (0,<f>). Therefore, when M is large, 
the ML estimator is computationally expensive, and requires to be solved numerically. 
Note that in contrast the ML estimator derived above does not involve any search 
procedure over the polarization vector. 

Eigenstructure-based techniques such as MUSIC and ESPRIT are computationally 
efficient for estimating source directions in a multi-source environment. These 
algorithms assume uncorrelated or partially correlated signals. In the following, a 
processing method is provided for utilizing eigenstructure-based techniques in the 
presence of fully correlated signals. 



Eigenstructure-based techniques for source localization, such as MUSIC, rely 
on identification of the signal and noise subspaces. In the presence of fully correlated 
signals the dimension of the signal subspace is smaller than the number of signals, M, 
and therefore, the signal subspace does not span the M-dimensional subspace of the 
spatial transfer functions, %{O X ^ X ^ X \...^{0 M ^ M ^ M ) . In this case, it is required to 
employ the information on the structure of the spatial transfer function, g(#,^,p), in 
order to determine the spatial transfer functions subspace. Spatial smoothing, forward 
backward averaging and redundancy averaging methods utilize the information on the 
structure of the spatial transfer function in order to estimate this subspace or part of it. 




The ML estimate of the signal can now be written as 




The Vector Sensor Averaging (VSS) Algorithm 
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The deficiency of the spatial smoothing method is the reduction of the effective array 
aperture length resulting in lower resolution and accuracy, while the deficiency of the 
redundancy averaging method is that its estimation errors bias does not vanish 
asymptotically for large numbers of measurements. In addition, both approaches are 
limited to the case of a linear equally spaced (LES) sensor array with far-field 
approximation. The forward backward averaging method assumes a symmetric array, 
far field approximation and unequal signal phases at the center of the array. 

In the methods of the present embodiments the vector sensor information is 
used in order to determine the subspace spanned by the steering vectors 
q(5„^)t...,q(fl w ,^), which enables estimation of direction of arrivals using eigen- 
structured-based methods, such as MUSIC. This objective can be obtained by Vector 
Sensor Averaging (VSS) method as described below. 

By substitution of eqn. (13) into eqn. (16) the measurement model at the array 
can be written in the form 

M 

y*=2h(^^)® A (^^)k^* +n *' k = i,...,K (3i) 

If we consider only the sensors of type / ( 1 </<£), . then the corresponding 
measurement vector, y lk , can be expressed as 

M 

y*=£fa(*-.rf.)® A i<*..*.)] ff-A+n* . k = l,...,K, / = 1,...,Z (32) 

m=\ 

which can be simplified to 

yft=Eq(^.^+»ftJ k = l,...,K l = l 9 ... 9 L (33) 

where A l {6 m ^ m ) is the /th row of the matrix A(0 m ,0 m ), z ml = A l (0 m9 ^ m )g m denotes 
the response of the /th type sensor for DO A (0 m ,<p m ), and n lk stands for the 
corresponding noise vector. 

Eq. (33) implies that each type of sensor array measurements provides a different 
linear combination of the vectors q(0 l9 ^ l \. ..,q(0 M ,<fi M ), as it would be the case for 
non-coherent signals. 

The information acquired by the L different sensor types helps to obtain a 
measurement space in which the signals are not fully correlated. We utilize this fact 
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in order to span the signal subspace, which is a necessary requirement of the eigen- 
structured based algorithms for source localization. 

Eq. (33) can be rewritten in matrix form as 

5 

y ft =Q(0,*>A+n ft , * = 1,...,A:,/ = 1,...,Z (34) 

where Q(0,&) = [q(0,MX-">4(0*' and z / =[ z u>- -> z miY - 
Therefore, the covariance matrix of each sensor type is given by 

R y = £[y tt yir] = ^ 2 Q^fQ"+R n , 1= i-L (35) 

10 in which cr s 2 denotes the signal power. The signal power is given by cr] - ^s k f\ 
and R nl denotes the corresponding noise covariance matrix. In this problem, the M 
signals are fully correlated. Accordingly, the signal covariance matrix of each sensor 
array type, cr* Qz t z?Q H is of rank one. In the VSS technique of the present 

embodiments, the covariance matrices , are smoothed for the L elements of 

15 the vector sensor. Consequently, the signal subspace is extended by averaging the L 
sensor type covariance matrices, i.e. 

1 L 

in which R z is defined asR 2 = —^x l z" . The rank of the new signal covariance 

matrix cr^QR 2 Q H is limited by min(rank(R z ),M) . Algorithms such as MUSIC, 
20 ESPRIT, etc, can use the corresponding sample covariance matrix, 

R = -— 2*-!^*^* 9 with stee " n g function q(#,^) for source localization. 

For determination of the M dimensional signal subspace, it is required that 
M < min(L, N) . This requirement can be alleviated if one can use other methods for 
signal decorrelation. For example, by applying the Forward-Backward averaging 
25 referred to above, the maximum number of the fully correlated signals, which can be 
localized, is doubled. 
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The vector sensor array contains NL sensors and therefore NL receivers are required 
for data collection. However the VSS computes the smoothed covariance matrix by 
averaging the NxN matrices {R y/ }^_ - In stationary scenarios the aforementioned 

matrices may be calculated for different periods, implying that one can use N 
receivers to collect the required data. Such may be achieved by use of a switching 
scheme. 

Reference is now made to Fig. 12, which is a simplified diagram illustrating 
switching between the elements. Switching control unit 200 controls switches 202 
provided at each of the sensor elements 1..4 to switch between the different dipole 
outlets at the signals. The signals detected as the switches cycle between the different 
dipoles are recorded sequentially into matrix positions ready for the subsequent 
processing stages. The use of switches in this manner greatly simplifies the input 
aspect of the system. 

Simulation 

Simulations were carried out to evaluate the performance of techniques according 
to the various above-described embodiments, that is to say for different vector sensor 
types. Performance was evaluated for different scenarios. Signal environments 
consisting of coherent multipath were simulated on the computer in order to verify the 
assertions of the previous sections. 

Simulations are presented for the following: 1) no preprocessing, 2) forward- 
backward averaging (FB), 3) VSS, and 4) VSS combined with FB (VSS-FB). The 
Cramer Rao bound (CRB) for source localization using a vector sensor array was 
derived in A. Nehorai and E. Paldi, "Vector-sensor array processing for 
electromagnetic source localization," IEEE Trans, on Signal Processing, vol. 42, pp. 
376-398, Feb. 1994, the contents of which are hereby incorporated by reference. 

The arrays used for simulations are twelve-elements linear arrays of vector 
sensors at half-wavelength inter-element spacing along the y-axis. Three kinds of 
arrays are used: 

a. Vertical polarized sensor array - Each antenna element consists of a 
vertical electric sensor (Fig. 3). This is the scalar sensor case and VSS 
preprocessing cannot be applied. 
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b. Dual polarized array - The vector sensors consist of vertical and 
horizontal electric dipoles (Fig. 4 with 5 = 0° according to fig.l 1). 

c. Quadrature polarized vector sensor array The vector sensors consist of 
four orthogonal components: two orthogonal electric dipoles and two 

5 orthogonal magnetic dipoles (sensors number 1, 2, 4 5 and 5 with 

S = 45° according to Eq. (11), assuming w = 1 ). 

d. In the simulations it was assumed that the sources were in the azimuth 
plane, that is 6 = 90° . 

In the first scenario, two equal power, fully correlated sources with DOA's 4°, 

10 0°, and elliptical polarizations Pl = (0.707e- j60 °, 0.707, p 2 = (0.707e j80 °, 0.707) were 
considered. The phase difference between the two incident signals at the origin was 
110° such that £ = (0.7.7e j50 ° 5 0.707e juo °), = (0.707e i80 °, 0.707). The number of 
samples taken from the array was 100. 

Reference is now made to Fig. 13, which is a graph that compares results of 

15 the above-described scenario for different decorrelation methods and also shows the 
Cramer Rao bound. The root mean square error (RSME) versus signal to noise ratio 
(SNR) is depicted. It can be seen that RSME of the FB-MUSIC decreases as the SNR 
increases, but it is not an efficient estimator even asymptotically. MUSIC with no 
preprocessing fails as expected with in the case of fully correlated sources. 

20 It is clear that VH MUSIC and FB MUSIC algorithms have the ability to 

decompose the coherent sources while the MUSIC algorithm by itself does not. The 
graph also shows that the ML can reach the Cramer Rao Bound (CRB). 

Reference is now made to Fig. 14, which is a graph presenting RSME versus 
the angle ?j 2 at SNR = 20dB for three different kinds of vector sensors, and three 

25 different types of preprocessing. The point A in Fig 14 is the point in Fig 13 at 
SNR = 2QdB . For this purpose the same scenario as used in Fig. 13 above is used with 
the SNR set to 20dB and the RMSE is plotted as a function of the polarization 
parameter rj 2 , see equation (6) above, where the subscript 2 refers to the number of 
the source. From Fig 14, it is clear that the FB MUSIC fails around two points, 

30 labeled B and C, which reflect cases in which the phase difference of the two signals 
at the array center is 0 or 180 degrees. At point D, where the polarizations of the two 
sources observed at the output of the array are proportional to each other, the VH 



29 

MUSIC fails. It is noted that in the above scenario, the absolute values of both 
polarization vectors were chosen to be identical. In practice, multi-path versions of 
the same signal do not usually have the same amplitude. In general cases, in which 
the polarization vectors are not identical, the limitation does not exist. The VSS FB 
5 MUSIC algorithm, which is a combination of VSS MUSIC and FB MUSIC, the 
above problem does not exist, since when the FB preprocessing fails to remove the 
singularity in the signal co variance matrix, it is removed by the VSS preprocessing 
and vice versa. The VSS-FB algorithm gives CRB resolution, as does the ML. 

Reference is now made to Fig. 15, which is a graph comparing four equal 

10 power, fully correlated signals with DOAs of -30, -10 0 and 20 degrees and randomly 
chosen polarization. The SNRs of all the signals were 15dB and 100 snapshots were 
collected from the array. The MUSIC and FB-MUSIC algorithms were applied to an 
array with vector sensors of type A, while the VSS and-MUSIC and VSS-FB-MUSIC 
algorithms were applied to an array of vector sensors of types B and C. As expected, 

15 one can observe that using an array of type A, MUSIC and FB-MUSIC fail to resolve 
the signal DOA's. sicne the VSS preprocessing in this case is unable to increase the 
rank of the signal covariance matrix to 4 as required to resolve four sources. By 
contrast, the VSS-FB-MUSIC with vector sensor of type B, and VSS-MUSIC and 
VSS-FB-MUSIC with vector sensors of type C are able to resolve the four fully 

20 correlated signals. 

Reference is now made to Fig. 16, which shows results using four additional 
coherent sources, that is eight sources, from directions 
-70 0 ,-50 0 ,-30 o ,-10 o ,0°,20 o 40 0 and 60°. From the figure we can see that only the 
VSS FB MUSIC with type c sensor, can resolve 8 coherent sources. All other cases 

25 fail. In the case of vector sensing using a quadrature array of type C, the VSS 
preprocessing results in a signal correlation matrix of rank 4. An additional FB 
preprocessing stage doubles the rank of the signal correlation matrix, in this case to 8, 
which is sufficient to allow the 8 sources to be resolved. 

It is appreciated that certain features of the invention, which are, for clarity, 

30 described in the context of separate embodiments, may also be provided in 
combination in a single embodiment. Conversely, various features of the invention, 
which are, for brevity, described in the context of a single embodiment, may also be 
provided separately or in any suitable subcombination. 
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Although the invention has been described in conjunction with specific 
embodiments thereof, it is evident that many alternatives, modifications and variations 
will be apparent to those skilled in the art. Accordingly, it is intended to embrace all 
such alternatives, modifications and variations that fall within the spirit and broad 
scope of the appended claims. All publications, patents and patent applications 
mentioned in this specification are herein incorporated in their entirety by reference 
into the specification, to the same extent as if each individual publication, patent or 
patent application was specifically and individually indicated to be incorporated herein 
by reference. In addition, citation or identification of any reference in this application 
shall not be construed as an admission that such reference is available as prior art to 
the present invention. 



